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I. INTRODUCTION 



Since the discovery of the equivalence between Af-fold super symmetry (SUSY) and weak 
quasi- solvability [1] in 2001, the research activity on exactly solvable Schrodinger operators 
has come into a new era. Indeed, we now recognize that all the exactly solvable ones 
whose solutions are expressible in terms of a classical orthogonal polynomial system are 
characterized as a particular case of type A TV-fold SUSY [2, 3] and that the ones which 
admit two linearly independent series solutions belong to type C A/"-fold SUSY [4]. In 
addition, a particular class of type B jV-fold SUSY constructed in Ref. [5] is now placed in a 
set of second-order linear differential operators preserving a so-called exceptional polynomial 
subspace of codimension one [6, 7]. A set of new quasi-solvable operators preserving an 
exceptional polynomial subspace of codimension two was also constructed in the framework 
of JV-fold SUSY [8], which is called type X 2 . For the development of jV-fold SUSY until the 
middle of 2000s and the terminology, see the review [9]. For its non-perturbative aspects, 
see also the review [10]. 

The mathematical structure of A/"-fold SUSY is characterized by a higher-order inter- 
twining relation if it is represented by differential operators. In general, however, it is quite 
difficult to solve directly the relation under consideration which actually consists of simul- 
taneous coupled nonlinear differential equations. Until now on, the most general form of 
A/"-fold SUSY is known only up to M = 4 [11]. They are all weakly quasi-solvable, that is, 
they keep a linear subspace annihilated by the corresponding A/"-fold supercharge invariant. 
It does not immediately mean, however, that some local solutions are available in a closed 
form. In other words, it does not automatically guarantee quasi-solvability in the strong 
sense. 

Hence, a natural question would arise; under what conditions an A/"-fold SUSY system 
admits a number of local analytic solutions. In the case where the number is two, it was 
proved [12] that a necessary and sufficient condition for it is to have type A 2-fold SUSY. 
The latter fact further enabled one to construct systematically shape-invariant potentials 
within this type of symmetry [13, 14]. These achievements clearly show the power of A/"-fold 
SUSY, and we expect that we would be able to step forward to more general cases, beginning 
with the case where the number of available analytic solutions is three. 

On the other hand, among the aforementioned four types of A/"-fold SUSY which are 
different with each other for M > 2, type A is the only one for which a necessary and 
sufficient condition is known. It seems that it is still quite hard to find it for an arbitrary 
Af > 2 even if we restrict our consideration to a particular type . Hence, it would be 
a realistic plan to make the examination from lower to higher values of Af step by step 
until when we eventually discover, if it exists, an inductive way toward a general statement 
applicable to an arbitrary Af > 2 case. 

In this paper, we examine under what conditions a one-body quantum Hamiltonian would 
admit three linearly independent local solutions in closed form. In addition, we investigate 
what type of 3-fold SUSY such a system must possess as a consequence of the equivalence 
between A/"-fold SUSY and weak quasi-solvability. We show that it is the most general type 
B, for which only a particular case has been explored. In other words, we prove that type B 
3-fold SUSY is a necessary and sufficient condition for a Hamiltonian to admit three linearly 
independent analytic local solutions. We find that there are eight linearly independent 
second-order linear differential operators of this kind. As a by-product of this finding, 
we obtain some new quasi-solvable operators preserving a three-dimensional monomial or 
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polynomial subspace which have not been considered in the literature. More precisely, we 
discover one new quasi-solvable operator for a type B monomial space, two for type C, and 
four for type X 2 . Furthermore, we find that all the eight quasi-solvable operators in the 
type A, type B, and type C cases are connected continuously via a one parameter. 

We organize the paper as follows. In the next section, we first construct the most general 
second-order linear differential operator which preserves a three-dimensional linear function 
space. In Section III, we solve directly the intertwining relation in type B 3-fold SUSY to 
obtain a necessary and sufficient condition for it. We then show that it leads to the operator 
which is identical with the one obtained in the previous section. The preceding two sections 
deal with particular cases of the most general type B SUSY. In Section IV, we consider 
the cases where the three-dimensional linear space is given by type A, type B, and type 
C monomial subspaces. In Section V, we first reformulate the most general type B 3-fold 
SUSY system in the previous section, and then apply it to the case of three-dimensional 
X 2 polynomial subspaces. In the last section, we summarize the results and discuss their 
implications and prospects for the future studies. 



II. QUASI-SOLVABLE OPERATORS PRESERVING A 3-DIM. SPACE 

Let us first consider a one-dimensional quantum Hamiltonian H which admits three 
linearly independent local analytic solutions in closed form to the corresponding Schrodinger 
equation 

HA(q) = E^q) (i = 1,2,3), (2.1) 

where each of does not necessarily normalizable. It is evident that the Hamiltonian H 
preserves a three-dimensional linear space V3 of functions defined by 

V3 = (Mq),Mq)Mq))- (2.2) 

Conversely, if a Hamiltonian preserves a linear space (2.2) spanned by three known functions, 
we can diagonalize it in the three-dimensional space algebraically and thus obtain three 
linearly independent local solutions to the corresponding Schrodinger equation. An available 
set of three functions differs depending on the Hamiltonian under consideration, but we can 
extract a common mathematical structure shared by all such systems. For this purpose, let 
us introduce a 'gauged' Hamiltonian H by H~ = tp^qjHipilq) -1 . Then, it immediately 
follows from the invariance HV3 C V3 that the gauged Hamiltonian H~ preserves the linear 
space 

V 3 "[*l = <W(*)>, (2-3) 
where the variable z and the function f(z) are defined by 

z = Mq)/Mq), /(*) = MqVMq)- (2-4) 

Here we note that it is not necessary for the function f(z) to admit an explicit form in terms 
of the variable z. Hence, any one-dimensional quantum Hamiltonian H in the physical q- 
space which has three linearly independent local analytic solutions can be converted with a 
gauge transformation and a change of variable z = z(q) to a gauged Hamiltonian H~ in the 
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gauged z-space which leaves the linear space (2.3) invariant. Such a gauged Hamiltonian H 
is no longer a Schrodinger operator but is a general linear second-order differential operator 
in the variable z: 

H- = -A(z)^-B(z)±-C(z), (2.5) 

where A(z), B(z), and C(z) are functions of z to be determined. To preserve the space 
(2.3), it is necessary and sufficient for H~ to satisfy 

H~\ = -C{z) = -c 2 f(z) - Cl z - c , (2.6) 
H~z = -B(z) - C{z)z = b 2 f(z) + b lZ + b , (2.7) 
H-f(z) = -A(z)f"(z) - B(z)f(z) - C(z)f(z) = a 2 f(z) + a,z + a , (2.8) 

where q, b iy and Oj (i = 0, 1, 2) are all constants. From these conditions, we obtain 

A(z)f"(z) = [(c 2 z + b 2 )f(z) + Cl z 2 + (6i + c )z + bo] f(z) 

- [c2f(z) + ciz + a 2 + co\f(z) - a x z - a , (2.9) 

B(z) = - (c 2 z + b 2 )f{z) - Cl z 2 - (co + h)z - b , (2.10) 

C(z) = c 2 f(z) + Cl z + c . (2.11) 

Substituting the latter formulas for A(z), B(z), and C{z) into (2.5), we see that H~ is 
expressed as 

H~ = - c 2 J 8 - ciJ 7 - b 2 J Q - (bi + c ) J5 - boJ^ 

+ (a 2 + c )J 3 + aiJ 2 + clqJi - c , (2.12) 



where J, (i = 1, . . . , 8) are given by 
Ji[z 



J 4 [z 
J 6 [z 

J 8 [z 



Id 2 z d 2 d 2 

J 2 [^]=^] = ^_" Mz]=f(z)J 1 [z]- J[ ~ } ' 



/"(^d^ 2 ' ZL J 11 J /"(^d^ 2 ' 7 U j ///( z ) dz 2' 

zf(z) d 2 d 

" /"(*) d^ 2 dz' 5 ^ " m ~ f"(z) d^ 2 *dz' 

ft w r 1 f( z )f'( z ) d2 // \ d ro 1^ 

= /(^N = - 7 ^ r ^-/(^, (2.13) 

zf\z)-f{z) e? 2 d 



and Jg is defined by 



T . . - /(*) d 2 d 1 , niA . 

j ' lz] = A»> d?-*di +i - (2 ' 14) 



Therefore, a linear space of non-trivial differential operators of at most second order which 
preserves the space (2.3) is spanned by the eight linearly independent operators Jj (i = 
1,...,8). 
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III. TYPE B 3-FOLD SUSY 



Type B A/"-fold SUSY was first discovered in Ref. [5] by a simple deformation of type A 
jV-fold supercharge. The component of type B 3-fold supercharge is given by 



) 



P 3 " = + W(q) - E(q) - F(q^j (J- + W(q)^ (J- + W(q) + P(g)) , (3.1 

where the three functions W(q), E(q), and F(q) are at present arbitrary. A pair of Hamil- 
tonians 
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H± =-\h +v±(q) < (3 - 2) 

is said to have type B 3-fold SUSY if it is intertwined by P 3 ~ in (3.1) as 

P?H* = H±If, (3.3) 

where P 3 + is the transposition of P 3 ~ in the g-space, P 3 + = (P 3 ~) T . A direct calculation shows 
that the intertwining relation (3.3) holds if and only if the potential terms in (3.2) have the 
following form 

V ± = -W 2 - - (2E' - E 2 ) - - (2F' + 2WF - 2EF -F 2 )±- (3W' - F') , (3.4) 
2 3 6 2 

and simultaneously the three functions W(q), E(q), and F(q) satisfy 

l q -E)F l - F - U - § ) = 0, (3.5) 



where Fi(q) and F 2 (q) are given by 



(3.6) 



F 1 = W' + EW - i (F' - 2WF + 2PP + F 2 ) , 
F 2 = E' + E 2 + ^ (F' - 2WF + 2EF + F 2 ) . 



(3.7) 



We note that the factorized form of the type B 3- fold supercharge component (3.1) is ex- 
panded as 



T = £ + (3-8) 

H k=0 H 



with 



w f = 3W-F, (3.9a) 

W W = 3W > + 2E' + 3W 2 - E 2 - 2WF - EF, (3.9b) 
W W = w" + E" + 3WW' + 2E'W - EE' - W'F - E'F 

+ W 3 - E 2 W - W 2 F - EWF. (3.9c) 
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Then, the same conditions (3.4)-(3.7) can be also derived from the general form of 3-fold 
SUSY in Ref. [11] by substituting (3.9) into the formulas in the latter reference. 

To solve the coupled differential equations (3.5) and (3.6), it is convenient to make a 
change of variable z = z(q) and to introduce a function f(z) defined by 

In terms of the new variable z and the function f(z), the two equations (3.5) and (3.6) are 
converted into 

^w-^(^w-nrl =0 - (3 - li; 
- 3 /'''(^^ w + ^^-m^ w _M |=0i ( , 12) 



dz 2f"(z)J 2W 2 V /"(*) f"(z) 2 J \ 1W 6 
where Fi(z(q)) = Fi(q) (i = 1,2). Eliminating F 2 (z) from (3.11) and (3.12), we obtain 



dz y">(z) 1 v ' f"{zf iv v v 1 /"'(*) 

This differential equation for Fi(z) can be integrated four times to yield 

F 1 (z) = d (zf(z) - 2f(z)) + C 2 f(z) + C 3 z + C 4 , (3.14) 

where Cj (i = 1,...,4) are integral constants. Substituting it into (3.11) and integrating 
the resultant differential equation for F 2 (z), we have 

F 2 (z) = -6(C lZ f'(z) + C 2 f(z) - C 3 z - C 5 ), (3.15) 

where C 5 is another integral constant. From the definition (3.10), we can express z and 
f(z(q)) in terms of E(q) and F(q) as 



;(g) = J dge/ 9d «'^'), /(*(,)) = J dq LfW) dq'e^' <*W')W)) 



(3.16) 



Substituting them into (3.14) and (3.15), we obtain expressions for F\(q) and ^2(9) defined 
in (3.7). Finally, substituting the resulting Fi(q) and F 2 (q) into (3.7), we obtain the entire 
relations to be satisfied among the three functions W(q), E(q), and F(q). For instance, the 
coupled differential equations (3.7) can be integrated for W(q) as 

W(q) = e-/ d "^) J dqer^'^') ^(q) + F M ~ E'{q) - E{qf ^ ^ 

Thus, we can express W(q) in terms of E(q) and F(q) by the substitution for the obtained 
Fi(q) and F 2 (q) into (3.17). Then, eliminating W(q) in (3.7) by using (3.17), we have a 
relation between E(q) and F(q). In this way, we can obtain the most general form of type 
B 3-fold SUSY. 
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In the next, we shall show that type B 3-fold SUSY systems are essentially identical with 
second-order linear differential operators preserving the space (2.3) investigated in Section II. 
For this purpose, let us first make a 'gauge' transformation on H ; 

H~ = c w *H~e- w z, (3.18) 

with the gauge potential W 3 (q) given by 

W 3 ± (g) = J dq(E(q) T W(q)). (3.19) 

Substituting (3.4) for V~(q) into (3.2), we have 

H- = -~ + {W{q) + E(q))± - Fi(q) - (3.20) 

where F 1 (q) and F 2 (q) are given by (3.7). Thus, if we make a change of variable z = z(q) 
and introduce two functions A(z) and Q(z) as 

2A(z(q)) = z\q)\ Q(z(q)) = -z'(q)W(q), (3.21) 
we can express the gauged Hamiltonian H~ in terms of z as 

Comparing it with (2.5), we immediately see the correspondence between them as 

B(z) = Q(z) - C(z) = h{z) + (3.23) 

On the other hand, we can convert the formulas among the functions of q in (3.7) into the 
following ones among the functions of z by employing (3.10) and (3.21) as 

p M n , (A 2r(z)A(z) + f"(z)(3A'(z)+2Q(z)) 

{ ) = ^ { ) ^ ' (3 24) 

p (A _ A „ (A , 2r(z)A(z) + r(z)(3A'(z) + 2Q(z)) [ " } 

F 2 (z) -A(z) + ^ . 

Substituting (3.24) into (3.14) and (3.15), we have coupled differential equations for A(z) 
and Q(z). The latter equations can be easily integrated out to yield 

2B(z) = 2Q(z) - A'{z) 

= AC lZ f{z) + 4C 2 /(z) - AC 3 z 2 - 2(C 4 + 3C 5 )z - 2C 6 , (3.25) 
A(z)f(z) = f'(z) [-2C lZ f(z) - 2C 2 f(z) + 2C 3 z 2 + (C 4 + 3C 5 )z + C 6 ] 

+ 2f(z)(C 1 f(z) - C 3 z - C 4 ) + + C 8 , (3.26) 

where Cj (i = 6,7,8) are additional integral constants. The explicit form of C(z) can be 
obtained by substituting (3.14) and (3.15) into the second formula in (3.23) as 

C{z) = -2CJ(z) + 2C 3 z + C 4 + C 5 . (3.27) 
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Comparing (3.25)-(3.27) with (2.9)-(2.11), we easily see that they are identical with each 
other and the parameters are related by 



c 2 = -2Ci, ci = 2C 3 , b 2 = -2C 2 , h = 2C 5 , b = C 6 , ^ ^ 

a,2 = C4 — C5, ai = — Cj, a = —Cs Co = C4 + C5. 

Therefore, the gauged type B 3-fold SUSY Hamiltonian (3.22) exactly coincides with the 
most general second-order linear differential operator (2.12) which leaves the linear space 
(2.3) invariant. The latter fact can be also derived directly by using the type B 3-fold super- 
charge without recourse to the explicit form of H~ . In the gauged z-space, the component 
of type B 3-fold supercharge reads as 

P 3 -W = e^-e-»- - £g) (3,9) 

and it actually annihilates all the elements of the linear space (2.3). In other words, 

V^=kerP 3 -. (3.30) 

On the other hand, it follows from the gauge-transformed version of the intertwining relation 
(3.3) with upper signs, namely, P%H~ = H + P%, that i/~kerP 3 ~ C kerP 3 ~. Hence, the 
equality (3.30) indeed means that any gauged type B 3-fold SUSY Hamiltonian H~ defined 
by (3.18) with (3.19) preserves the linear space (2.3). 

One of the most advantageous and powerful aspects of the framework of A/"-fold SUSY is 
that we can construct simultaneously another weakly quasi-solvable Hamiltonian H + which 
is almost isospectral with H . According to Ref. [4], the partner Hamiltonian H + of 3- 
fold SUSY in the physical g-space is connected to the gauged z-space by another gauge 
transformation induced by the gauge potential W^(q) already defined in (3.19). With the 
latter gauge transformation, H + is transformed to a second-order linear differential operator 
H + of the variable z which must have the following expression: 

H + = e w 3 + tf+ e -w 3 + 

= - A{z)— 2 + ^-ii + Q(z)j - - C{z) - 2Q\z) 

+ A\z)wf\z) + 2A(z)wf ] '(z), (3.31) 

where A(z) and Q(z) are the same functions as the ones in (3.21), C(z) is the zeroth-order 
term of H~ in (2.5) and thus is given by the second formula in (3.23) in our type B case, 
and w 2 3 \z) is defined via the expanded form of the gauged 3-fold supercharge component 
P 3 ~ under consideration: 

V fc=0 / 

Thus, in our type B case (3.29), it is given by 

4 3, M = -fjf- (3.33) 
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Substituting (3.25)-(3.28) and (3.33) into (3.31), we can express H + as 

H+ = - c 2 K s - Cl K 7 - b 2 K 6 - (61 + c )K 5 - b K 4 

+ (a 2 + c )K 3 + a x K 2 + a K 1 - c , (3.34) 

where the operators Ki (i — 1, ... ,8) are given by 



K r 1 1 d 2 d - r\ z f 

lN /" (2) d^ 2 f"{zf dz f"{zf 
K 2 [z] = zKM\ - KMl KM] = f{z)KM\ ~ f'(z)K [z] + 1, 

KM] = f{z)KMl KM] = f'(z)K 2 [z], KM] = f'(z)KM], 

KM] = (zf'(z) - f(z))KM\, KM] = (zf'(z) - f(z))K 3 [z], 

and K Q is defined by 



(3.35) 



K ^ = m^ + 7W (3 - 36) 



By applying the formulas (3.7), (3.10), (3.21), (3.23), and (3.33), the gauged partner Hamil- 
tonian H + in (3.31) is expressed in terms of q as 

H + = - + (E- W)^- + W + EW 
2 dq 2 dq 

- - (E' + E 2 + 5F' + 2WF - 2EF - F 2 ) . (3.37) 
6 

Then, we can easily check that the potential term in the original Hamiltonian H + = 
e~ w s H + e w 'i coincides exactly V + (q) in (3.4), which was derived from the direct calculation 
of the intertwining relation (3.3). 

The solvable sector V 3 preserved by H + is characterized by the gauge-transformed oper- 
ator of the transposed component P 3 + of type B 3-fold supercharge 

*=^-* = -*<Hs + ?w)- (3 - 38) 

In fact, the gauge transformation with W 3 of the intertwining relation (3.3) with lower signs, 
namely, P 3 H + = H~P 3 + , immediately leads to H + ker P 3 + C kerP 3 + . Hence, we obtain from 
(3.38) 

Vt\z] = kerP+ = ^(lJMlzf'iz) - f(z)). (3.39) 

We can easily check that the eight linearly independent operators K^ {% — 1, . . . , 8) in (3.35) 
indeed preserve the linear space (3.39). It is worth noting that the space V3" in (3.39) is 
related to the one V 3 in (2.3) by 

Vi\z] = f\z)- l %[w]\ w=nz)J{w)=zfl{z) _ f{z) . (3.40) 
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Then, it follows from this relation that a set of eight linearly independent second-order 
linear differential operators preserving can be constructed directly from Jj (i — 1, ... ,8) 
in (2.13) as 

/"(,)- 1 J l Hf(,)|^, (z))/w ^ w _ /w (i = 1,...,8). (3.41) 

In fact, we can check that all the operators K^ {% — 1, . . . , 8) in (3.35) are expressed as linear 
combinations of the ones in (3.41) as 

km = rw^JMnz), km] = nzr'Mwwxz), 

KM] = f"(z)-\j 5 [w] - J 3 [w] + !)/"(*), KM] = fXzr'JMU'Xz), , 42) 

km] = rw^JMn*), km] = nzy^wiz), [ ■ j 

tf 7 [z] = /"(^-VeH/'^), km = nzMMwW'iz), 

where the substitution w = f'(z) and f(w) = zf'(z) — f(z) in each operator Ji[w] has been 
understood. 



IV. TYPE A, B, AND C 3-FOLD SUSY WITH MONOMIAL SUBSPACES 

In this section, we shall examine some particular cases of the most general type B 3-fold 
SUSY investigated in the previous section. Virtually all the A/"-fold SUSY models (including 
the ones called with other terminologies) so far constructed in the literature are of the types 
which preserve a monomial or polynomial subspace. The latter types of second-order linear 
differential operators for the monomial cases were first classified in Ref. [15]. According to 
it and to the later reconsiderations in Refs. [4, 5], there are essentially three inequivalent 
monomial subspaces, except for a few low-dimensional ones, preserved by second-order linear 
differential operators. They are respectively called type A, type B, and type C monomial 
subspaces [4, 9]. A/"-fold SUSY systems associated with them were investigated in general 
fashions in Refs. [2, 3, 16] for type A, in Ref. [5] for type B, and in Ref. [4] for type C 
monomial subspaces. 

As we have shown in the previous section, any second-order linear differential operator 
which preserves a three-dimensional linear space of functions admitting an analytic expres- 
sion has type B 3-fold SUSY. It means that all the aforementioned systems which preserve a 
monomial or polynomial subspace must have type B 3-fold SUSY too when Af = 3. Hence, 
it is interesting to know how those systems in the literature would be realized in the most 
general type B 3-fold SUSY obtained in the previous section. 

Let us first begin with type C since, as we will see later, type A and type B can be 
regarded as special cases of type C. A type C monomial subspace Vjy/- ^ of dimension M is 
defined by [4] 

V'SU = <!> *, • • • , z\ z X+ \ ■ ■ ■ , z X+N2 Mi A/" = A/"i + A/" 2 > 3, (4.1) 

where M\ and A/2 are positive integers and A is a real number. The linear space (2.3) we 
have considered reduces to a type C monomial subspace of dimension 3 with J\f\ = 2 and 
M 2 = 1 when f(z) = z x : 

V s -| /W ^ = (1, *, ^ A > = Vg } (A ^ -2, -1, . . . , 3). (4.2) 
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The corresponding gauged 3-fold supercharge component (3.29) which annihilates this space 
(4.2) reads as 

fi n ^( d2 A-2\ d 2 _ f3l/ , A - 2 



dz 2 z / dz 2 ' 

where iu 2 ^(;z) is calculated with (3.33). The set of eight second-order linear differential 
operators (2.13) which preserve the three-dimensional type C monomial subspace (4.2) is 

d^ 2 ' 

(4.4) 



A(A 


- 1) Ji : = 


= Ji(A) = ^" A ^, A(A - 1) J 2 := J 2 (A) = 


A(A 


- 1) Js ■-- 


= J 3 (X) = z 2 ^ (A - 1) J 4 := J 4 (A) = ^ 


(A- 


1)J 5 : = 


MX) = z ^-(x-i )z ± 


(A- 


1) J 6 := 


Ux) = z ^-(x-iyl 


XJ 7 


:= Jr(A) 


= z 3 — — - — Xz 2 — — h Xz, 
az z dz 


AJ 8 


:= J 8 (A) 


d 2 d 

az 2 az 



On the other hand, in the existing literature there are six linearly independent linear differen- 
tial operators of at most second order preserving the linear space (4.2), which are summarized 
in (Al) and (A2). The correspondence between (4.4) and the latter is as follows: 

(A-1)J( C) (A) = J 3 (\)- MX), J (C) (A) = J 4 (A), J( C) (A) = MX), 
jf } (A) = MX), jfl(x) = MX), J$(\) = MX). 

Hence, we see that the two operators Ji(A) and J 2 (A) have been missed so far in the literature. 
By using (4.4) and (4.5), the most general second-order linear differential operator (2.12) 
having an invariant type C monomial subspace (4.2) is expressed as 

H- = - 52J s( A) - ^(A) + ^A(A) + j^MX, 

- Sjjfo'tA) - S 2 J ( °\\) - ii J<?(A) - h 4 c> (\) - co, (4.6) 

where 

° 3= A' ° 2 = T^T - A(A^T) ' ai= A^T' " 6 o = ^ + c o- ^ 

When we put c 2 = b 2 = ai = a = 0, it reduces to the gauged type C 3-fold SUSY 
Hamiltonian constructed in Ref. [4], Eq. (3.8). 

The gauged 3-fold superchage component (3.38) and the solvable sector (3.39) annihilated 
by it for the plus component in this case read as 

* = + W = *<M,^>. (4.8) 
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(4.9) 



The set of eight linearly independent second-order linear differential operators (3.35) which 
preserve the latter space reduces to 

A(A - 1)*! := K^X) = z 2 ~ x ^- 2 + (A - 2)^~ A ^ - (A - 2)z~\ 

A(A - 1)K 2 := K 2 {X) = z 3 ~ x -^ + (A - 3)^" A ^ - 2(A - 2)z 1 ~\ 

A(A - 1)K S := Ks(\) = z 2 ^ - 2z± + 2, 

(A - l)K A := K 4 (X) = + (A - 2)^- - (A - 2)z-\ 

(A - 1)K 5 := K 5 (X) = z 2 -^ + (A - 3)^ - 2(A - 2), 

(A - 1)K 6 := K 6 (X) = z x+1 -^ - 2z x ^- + 2z x -\ 

dz z dz 

\K 7 := k 7 (X) = z^ 2 + (A - 3)z 2 ± - 2(A - 2)z, 

XK 8 := K 8 (X) = z x+2 ^ - 2z x+1 ± + 2z x . 

dz 1 dz 

The relation to the known quasi-solvable operators listed in (A3) is as follows: 

(A - 1)^ C) (A) = K 5 (X) - K 3 (X) + A - 1, K^(X) = K 4 (X), 
K^(X) = K 3 (X), K < S(X) = K 7 (X), Kf}(X)=K 1 (X), K$(X)=K 2 (X). 



(4.10) 



Hence, we see again that the two operators Kq(X) and K 8 (X) have been missed in the 
literature. 

Finally, we note that the function F(q) defined in (3.10) is calculated in this case as 

and coincides (up to a multiplicative factor) with the ad hoc constraint F(q) = z'(q)/z(q) 
made in Ref. [4], Eq. (3.17). We now understand that it naturally comes from the general 
relation (3.10) for the general case of an arbitrary f(z). 

Let us next consider the case associated with a type B monomial subspace. An M- 
dimensional type B monomial space is defined by [4, 5] 

V#> = <M,..., ^" 2 ,^>. (4.12) 

For M = 3, it can be realized as a special case of the linear space (2.3) with f(z) = z 3 or as 
a special case of the type C monomial space (4.2) with A = 3: 

%\ m ^=^\^ 3 = (l,z,z 3 ) = vP. ( 4 - 13 ) 

The corresponding gauged 3-fold supercharge component (3.29) which annihilates this space 
(4.13) reads as 

* = 1) * H w = 4 (414) 
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where w^\z) is calculated with (3.33). The set of eight second-order linear differential op- 
erators (2.13) which leave the three-dimensional type B monomial subspace (4.13) invariant 
is 

j.<3> = 4(3) = £, A(3)=^, 

aw - 4 - 4- J »< 3 > - ^ - H? J «< 3 > - ^ - 2 4- (4 - l5) 

J 7 (3) = ^ - 3^iL + 3*, J„(3) = - + 32 3 . 

dz 1 dz dz A dz 

On the other hand, in the existing literature there are seven linearly independent linear 
differential operators of at most second order preserving the linear space (4.13), which are 
summarized in (Bl) and (B2). The correspondence between (4.15) and the latter is as 
follows: 

2jf) = J 3 (3)-J 5 (3), ji B) = J 2 (3), jf = J 4 (3), 
jff = J 3 (3), jf } = J 7 (3), j£! = J«(3), jg} = J s (3). 

Hence, we see that the one operator Ji(3) has been missed so far in the literature. By using 
(4.15) and (4.16), the most general second-order linear differential operator (2.12) having 
an invariant type B monomial subspace (4.13) is expressed as 

ff- = -|j 8 (3) + ^Ji(3)- aijJ^-boJ^-co, (4.17) 

i,3=+,0,-, 
i>3 



where 



2a ++ = b 2 , 3a +0 = Ci, 6a o = — ^2 + 36i — 2c , 
2ao_ = &o, 6a = — ai, — 6q = &i + cq. 



(4.18) 



When we put c 2 = ao = 0, it reduces to the gauged type B 3-fold SUSY Hamiltonian in 
Ref. [5], Eq. (3.12). 

The gauged 3-fold superchage component (3.38) and the solvable sector (3.39) annihilated 
by it for the plus component in this case read as 



(4.19) 



The set of eight linearly independent second-order linear differential operators (3.35) which 
preserve the latter space reduces to 



*a(3) = ^-2z< + 2, /v ; (3) : 

X 5 (3) = z^ + 2, tf 6 (3) = /-^ - 2, 3 A + 2z\ 

K 7 (3) = z 3 -^ - 2z, K 8 (3) = z 5 ^ - 2/A + 2 z\ 
dz 1 dz 1 dz 



(4.20) 
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The relation to the known quasi-solvable operators listed in (B3) and (B4) is as follows: 

= K 5 (3) - 2, tfg> = jf 7 (3), K™ = K 9 (3), K ( ^=k 8 (3). 

Hence, we see again that the one operator Ki(3) has been missed in the literature. 

Finally, the function F(q) defined in (3.10) is calculated in this case as F(q) = z'(q)/z(q) 
and coincides exactly with the ad hoc constraint on F(q) made in Ref. [5], Eq. (3.6). As in 
the case of type C, it naturally comes from the general relation (3.10). 

In the last, we shall consider the case associated with a type A monomial subspace. An 
M- dimensional type A monomial space is defined by [2, 4] 

V#> = <1, a"" 1 ). (4.22) 

For M = 3, it can be realized as a special case of the linear space (2.3) with f(z) = z 2 or as 
a special case of the type C monomial space (4.2) with A = 2: 

Vs1 /W= * = V£?| A=2 = < W) = V? } - (4.23) 

The corresponding gauged 3-fold supercharge component (3.29) which annihilates this space 
(4.23) reads as 

P 3 - = ^(g) 3 -^, 4 31 (*) = 0, (4.24) 

where wf\z) is calculated with (3.33). The set of eight second-order linear differential 
operators (2.13) which preserve the three-dimensional type A monomial subspace (4.23) is 

4- /•<'>=*■;& -4- (4 - 25) 

J 7 (2) = ^ - 2z* A + 2 2 , J 8 (2) = 2 4 Al - 22 s A + 2 J 2 . 
dz 2 dz dz 2 dz 

In contrast with the previous two cases of type C and type B, it is entirely equivalent with 
the set of eight linearly independent linear differential operators of at most second order 
preserving the linear space (4.23) already appeared in the existing literature, which are 
summarized in (CI). The correspondence between (4.25) and (CI) is as follows: 

J™ = J 2 (2) - J 4 (2), 4 A) = J 3 (2) - J 5 (2), Jf) = J 6 (2) - J 7 (2), 
A A } = J l (2), ^ A) = J 2 (2), 4 A) = M3), = J 6 (3), ji A + } = J 8 (3). 

By using (4.25) and (4.26), the most general second-order linear differential operator (2.12) 
having an invariant type A monomial subspace (4.23) is expressed as 

H~ = - E W A) "Cb. ( 4 - 2 7) 

«j'=+,o,- »=+,o,- 
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where 

2a ++ = c 2 , 2a +0 = 2b 2 + c 1 , 2a 00 = -a 2 + 26i + c , 

(4.28) 

2a _ = — ai + 2b , 2a = — a , 26 + = ci, 6 = h + c , o_ = bo- 
lt exactly coincides with the gauged type A 3-fold Hamiltonian in Ref. [2], Eq. (4.22), which 
is originated from the sl(2) Lie-algebraic quasi-solvable operators in Ref. [17]. 

The gauged 3-fold superchage component (3.38) and the solvable sector (3.39) annihilated 
by it for the plus component in this case read as 

A + = -^(<?) 3 ^, V 3 + = <M,* 2 >, (4.29) 

and are entirely identical with the ones for the minus component, (4.23) and (4.24). The set 
of eight linearly independent second-order linear differential operators (3.35) which preserve 
the latter space reduces to 

* I( 2) = £, &<2) = *,(2) = ^*-24 + 2, 

= A- 5( 2)^- 2 il, (4.30) 

and is again equivalent with the set of eight linearly independent linear differential operators 
of at most second order preserving the linear space (4.23) in (CI). The correspondence 
between (4.25) and (CI) is as follows (note that k\ A) = j[ A) and K (A) = j\ A) = +, 0, -) 
in type A): 

K [ A) = K 4 (2) - K 2 (2), K iA) = K 5 (2) - K 3 (2) + 2, 

K {A) = K 7 (2)-K 6 (2), K i _ A l = K 1 (2), K {A) = K 4 (2), (4.31) 
K$> = 2K 5 (2) - K 3 (2) + 2, K$> = K 7 (2), K$ = K 8 (2). 

Finally, the function F(q) defined in (3.10) vanishes in this case, F(q) = 0, which naturally 
explains the reason why type A 3-fold SUSY can be characterized only by two functions 
E(q) and W(q). 



V. TYPE X 2 3-FOLD SUSY WITH POLYNOMIAL SUBSPACES 



In a certain problem, it is more convenient to start with a more general three-dimensional 
linear space of functions in a gauged space than (2.3). Type X 2 3-fold SUSY [8] is such an 
example. Hence, we shall first reformulate type B 3-fold SUSY in a more general setting. 
Instead of the linear space (2.3), we shall consider as an invariant subspace 

Vs~M = (<Pi{u),<P2(u),<P3(u)), (5-1) 
in a gauged w-space. The relation to the linear space (2.3) in a gauged z-space is obviously 

V 3 -[u] = ^i(«)V 3 -N| z= - 2(u) ^ i(u))/(z)= - 3(u) ^ i(u) . (5.2) 
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With this relation, we can easily produce the most general type B 3-fold SUSY in the -u-space 
associated with the linear space (5.1). Indeed, it follows from (5.2) that the gauged type B 
3-fold supercharge component P^~ [u] and the eight linearly independent second-order linear 
differential operators J'i[u] (i — 1, ... ,8) which preserve the linear space (5.1) in the gauged 
«-space are constructed from the ones in the z-space, namely, P^[z] in (3.29) and J^z] in 
(2.13) via the relations 



\Z=<p2(u)/ipi(uy 



P'f[u] = (p 1 {u)P^[z](p 1 {u)- 1 \ 
Their explicit forms are respectively given by 



(5.3) 



N - u (q) ^ + ^ + mi mi2 J ^ + _ i W J ^ „ J , (5.4) 

where each of the function W^j{u) (Wij t ki(u), respectively) is the Wronskian of <Pi(u) and 
<fij(u) (of Wij(u) and W k> i(u), respectively) 

Wij(u) = <j?i(u)(pj{u) - (pi(u)<fj(u), (5.5) 
W ijM (u) = W'^uWUu) - W t ,(u)W^(u), (5.6) 

and by 

, M _ ( d 2 d w^r - w 2 ^'[ 

J ll u \ - — 77} 17?2 ~ T77--T-. + 



W 3 i,2i \du 2 W 2 ,idu W 2 ,i<fii 

J 4l U \ — TT7 I JO MT J„. + 



(5.7) 



^31,21 \du 2 W 2:1 du W 2j i(f>! J W 2: i \du fix 

J'M] = ^J[[u], J'M = ^J[[u], J'M] = ^J'M, 
V?i Vi fi 

J'M] = ^J'M, J'M = ^*], J'M] = ^J'Ml 
<pi <pi <pi 

where J'M\ ls defined by 

, , , = W^f f^_Wi 1 d_ W^&-W 2 tf{ \ _^fi(d__&\ 

9U W 31 ,2i \du 2 W 2A du W 2 ,i<pi J W 2A \du <pj { ' } 

The plus component of the gauged type B 3-fold supercharge -Ps + [w] and the linear space 
V3 + [w] annihilated by it in the -u-space are related to the ones in the z-space, namely, to 
P^[z] in (3.38) and to V£[z\ in (3.39) respectively as 



P^[U] = ^ 1 ( M ) 3 ^ 2 , 1 ( M )- 2 P3 + [^1U 2 , 1 ( M ) 2 ^ 1 H- 

v^M = ^iH 3 ^,iH- 2 V3 + W|^ 2(u)M(u) - 

Their explicit forms are respectively given by 



\z=lf2(u)/lf>l(u) ' 



(5.9) 



P>+\u] - —u'(q) 3 { ± + &\(±_& + 3A(±_&^ + ^) (5 10) 
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and by 



v 3 + M = 



W 3 i,2l(«) 



(W 2A (u) : W 3A (u),W 3 , 2 (u)). 



(5.11) 



The eight linearly independent second-order linear differential operators K'^u] (i — 1, . . . , 8) 
which preserve the latter space (5.11) are also constructed from Ki[z] in (3.35) via the 
relation 



K>[u] = ^ 1 ( M )^ 2 , 1 ( M )- 2 ^[^^2,i(«)^i(«)- 3 L = ^ (tt)M(tt) 
and are explicitly given by 

d 2 /2# W^ 21 \ d <p\ 



(5.12) 



K[[u] 



^31,21 
^31 21 



w- 



31,21 



d-u 2 



^31,21 / d-u v^i W"2,i 

_ W^l + ^31,21 

£l ^2,1 ^31,21/ ^31,21. 

<p3 r W r ! W3,l . 



f - ^[«], = f - ^k'M + 1, 

Vl Vl ^2,1 



(5.13) 



WW 



W 2) i" 



^2,1 



3 ["J) 



VV 2 l VV 2 i 



where K' Q is defined by 



K >\ 1 _ (^2,i) 2 / d ft W' %x W^ 21 

° m w 31>21 \du q> x w 2 ,i ^31,21 



(5.14) 



With these preliminaries, let us construct the most general type X 2 3-fold SUSY system. 
The three-dimensional type X 2 polynomial subspace is introduced as (5.1) with (p n {u) (n = 
1, 2, 3) given by [8] 



(p n {u] a) — (a + n — 2)u n+l + 2(a + n — l)(a— l)u n + (a + n)(a — l)au 



n-l 



(5.15) 



where a is a free parameter. The gauged type X 2 3-fold supercharge component is calculated 
by using the formula (5.4) as 



(- 


f'a+3 ' 




( — 


fa+2 i 


f r:( 




fa+3j 


Ja+1 









P^[u-a}=u'{qfj o 

where the function f a {u) is given by 

fa( u ) — f( u 'i a) = u 2 + 2(a — l)u + (a — l)a. 



la+l 



(5.16) 



(5.17) 



It exactly coincides with Eq. (3.6) for Af = 3 in Ref. [8]. However, the latter reference 
reported only four linearly independent second-order linear differential operators j\ X2 \a) 
(i — 1, . . . ,4), summarized in (Dl), which left the type X 2 polynomial subspace (5.1) with 
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(5.15) invariant. A direct calculation of the eight operators J-[w] in (5.7) using (5.15) shows 
that the former are expressible as linear combinations of the latter as 

j[ X2 \a) = 2(a + 3) J 2 [u; a] - 2J' 3 [u; a] - (a + 2)J' A [u; a] + J' 5 [u; a] - 2, 
a(a + 1) J {X2 \a) = 2a{a + 3) (a + 2){a - l)J[[u; a] 

- 2a{a - l)(3a 2 + 12a + 13)J 2 [u; a] + 2(a 2 - 3a - 2)J' 3 [u; a] 
+ (a + 2)(a - l)(3a 2 + 6a + 4)J' 4 [u; a] + 4(a + 2)J' 5 [u; a] 

- a 2 Jg[w; a] + 2(a 2 - l)J 7 [u; a] + 2(a - 2)(a 2 - 1), 

J (X2 \a) = 2(3a 2 + 7a + 6) J' 3 [u; a] - (3a 2 + 5a + A)J' b \u; a] ^ ^ 

+ aJ^u; a] - 2(a - l)J' 7 [u\ a] + 4(a + 4)(a - 1), 

2(a + 1) j[ X2 \a) = - 4a(a + 3) 2 (a 2 - l)J' 2 [u; a] - 2(a - l)(7a 3 + 31a 2 

+ 54a + 36) J 3 [u; a] + 2a(a + 3)(a + 2) 2 (a - l)J' 4 [u; a] 
+ (a - l)(7a 3 + 31a 2 + 54a + 48) J' 5 [u; a] - (3a 3 - 5a 2 

- 14a - 8) Jq[u; a] + 2(a - l)(3a 2 + a - l2)J' 7 [u- a] 
+ 4(a - 1) Jg[u; a] - 8(a - l)(a 3 + 7a 2 - 10). 

Hence, there are four other linear combinations of J-[u; a] (i = 1, . . . , 8) which are linearly 
independent of J^ X2 \a) (i = 1, . . . , 4). 

The plus component of the gauged type X 2 3-fold supercharge P 3 + is similarly calculated 
from the formula (5.10) with (5.15) as 

p, + [u] = u , { f fd + fjx+A l^±f± + !f±l \ f -f^(^ + (5.19) 

which exactly coincides with Eq. (3.18) for J\f = 3 in Ref. [8]. The linear space V3 + [w] (5.11) 
annihilated by it reads by using (5.15) as 

V 3 + [u; a + 3] = f a (uY l f a+z {uy l {xi{u, a + 3), x 2 («, a + 3), x 3 (w, a + 3)) 

= f a {u)~ 1 faM- 1 % X2b) [u-a + ^ (5.20) 

where the function x n (w; a) (n — 1, 2, 3) is defined by 

Xn(u; a) — (a - n)(a - n + l)u n+1 + 2(a - n - l)(a - n + l)(a - l)u n 

+ {a-n-l)(a-n)(a-l)au n -\ (5.21) 

and thus reproduces Eq. (3.20) for Af = 3 in Ref. [8]. It is evident that there are eight 
linearly independent second-order linear differential operators K-[u;a + 3] which preserve 
the polynomial subspace V 3 X2b \u; a + 3] appeared in (5.20). Since K'^u] (i = 1, . . . , 8) in 
(5.13) leave the linear space V' 3 + [u; a + 3] invariant, it follows from (5.20) that they are given 
by 

k[[u- a + 3] = f a {u)- l f a+ ^r l K[[u- a + 3]/ a+3 («)/ a («) (i = 1, . . . , 8). (5.22) 

There were again only four linearly independent second-order linear differential operators 
K\ X2 \a) in Ref. [8], summarized in (D2), which preserved the linear space V^ 2h \u\a}. A 
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direct calculation of the eight operators K[[u] a + 3] in (5.22) using (5.13) and (5.15) shows 
that the former are expressible as linear combinations of the latter as 

K {X2 \a) = 2aK' 2 [u- a] - 2K' 3 [u; a] - (a - l)K' 4 [u; a] + K' 5 [u; a] - 2, 
(a - 2) (a - 3)K {X2 \a) = 2a(a - l)(a - 3) (a - 4)tf([u; a] - 2(a - 3) (a - 4) 

x (3a 2 - 6a + 4)K' 2 [u; a] + 2{a 2 - 9a + 16)K' 3 [u; a] 
+ (a - l)(a - 4) (3a 2 - 12a + 13)K' 4 [u; a] 
+ 4(a - l)K' 5 [u; a] - (a - 3) 2 A^[w; a] 
+ 2(a - 2)(a - 4)K' 7 [u; a] + 2(a - 2) (a - 4)(a - 5), 

K {X2 \a) = 2(3a 2 - 11a + 12)K' 3 [u; a] - (3a 2 - 13a + 16)K' 5 [u; a] 

+ (a - 3)K' 6 [u; a] - 2(a - 4)K' 7 [u; a] + 4(a + l)(a - 4), 
2(a - 2)K[ X2 \a) = - 4a 2 (a - 2) (a - 3) (a - A)K' 2 [u; a] - 2(a - 4) (7a 3 - 32a 2 
+ 57a - 3Q)K' 3 [u; a] + 2a(a - l) 2 (a - 3)(a - ^K'^u; a] 
+ (a - 4) (7a 3 - 32a 2 + 57a - 24)K' 5 [u; a] - (3a 3 - 32a 2 
+ 97a - 92)K' 6 [u; a] + 2(a - 4) (3a 2 - 17a + 12)K' 7 [u; a] 
+ 4(a - 4:)k' 8 [u; a] - 8(a - 4)(a 3 - 2a 2 - 15a + 26). 

(5.23) 

Hence, there are four other linear combinations of K[[u\ a] (i — 1, ... ,8) which are linearly 
independent of K\ X2 \a) (i = 1, . . . , 4). 

VI. DISCUSSION AND SUMMARY 

In this paper, we have shown that type B 3-fold SUSY is a necessary and sufficient 
condition for a Schrodinger operator to admit three linearly independent local analytic 
solutions and thus to preserve a three-dimensional non-polynomial linear space. The most 
general operator of this kind consists of eight linearly independent differential operators of 
at most second order. As a by-product of the latter finding, we have found that all the 
known 3-fold SUSY associated with monomial or polynomial subspaces, which are regarded 
as particular cases of the most general type B 3-fold SUSY, also admit eight operators 
for each component Hamiltonian, some of which have been missed in the literature. In 
particular, the three different types of 3-fold SUSY associated with monomial subspaces, 
namely, type A, type B, and type C, are connected continuously via a one parameter. The 
results obtained here are full of implications for the future studies, and we would like to 
close this paper by referring to some of them. 

The fact that the number of linearly independent quasi-solvable operators of at most 
second order preserving a three-dimensional linear function space is always the same casts 
a natural question: is it true for any dimensional linear function space? According to the 
classification of quasi-solvable operators which leave a monomial space invariant done in 
Ref. [15], there are no extra operators in the cases of type B and C monomial subspaces 
(the general case (4) and the special case C and D there) except for the three-dimensional 
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and a particular four-dimensional ones. The number of linearly independent quasi-solvable 
operators of at most second order preserving a type A monomial subspace (the general 
case (3) in Ref. [15]), on the other hand, is always eight regardless of the dimension of the 
subspace. Hence, at the first sight, it does not seem to be true in the cases of function spaces 
of more than three dimension. 

Let us next observe the problem from a different point of view. A crucial consequence of 
the fact in the present three-dimensional case is that the eight linearly independent quasi- 
solvable operators preserving monomial spaces of type A, type B, and type C are connected 
continuously with each other via a one parameter. In this respect, we note that these three 
different types are connected continuously, at least at the level of monomial subspaces, also 
in any dimension of more than three. In fact, an A/"-dimensional type C monomial subspace 
(4.1) with Mi = Af - 1 and Af 2 = 1 reads as 

V^ M = (1,,,...,^- 2 ,A (6.1) 

Evidently, it reduces to an Af- dimensional type B monomial subspace when A = Af and to 
a type A one when A = Af — 1. Hence, in contrast with the previous observation, it seems 
that the latter two types are connected continuously with the particular case of type C via 
the parameter A also in any dimension Af > 3. 

If it is indeed the case, the number of linearly independent quasi-solvable operators in the 
cases of type B and type C must be the same as the one in the type A case, namely, eight. 
If it does not the case, on the other hand, some irregular phenomena would take place in 
the limits A — > 1 and 2. It could be superficial, however, arising from the restriction to at 
most second-order linear differential operators. Actually, if we restrict our consideration to 
first-order ones in the present Af = 3 case, they are no longer connected continuously with 
each other and the number of linearly independent operators varies (one in type B and type 
C, three in type A). Hence, it could be the case that the number is constant for operators 
of a certain order higher than two in a linear function space of dimension higher than three. 

Another intriguing fact is that there is at least one extra quasi-solvable second-order 
operator for a particular four- dimensional type C monomial subspace with Af\ = Af 2 = 2 
(the special case B in Ref. [15]): 

V^ = (l,z,z\z x+1 ). (6.2) 

Hence, we expect that a detailed examination for the four-dimensional case would give us 
further insights into the issue. 

The fact that the most general type B 3-fold SUSY system admits three linearly indepen- 
dent local solutions of non-polynomial type in closed form indicates that there would be a 
possibly tremendous number of Schrodinger operators having non-polynomial local solutions 
which have not been found yet. Until now, virtually all the Schrodinger operators known 
to have local analytic solutions and thus be (quasi-)solvable are those whose solutions are 
expressible as products of a polynomial in a particular function and a gauge factor (the 
ground-state wave function when exactly solvable). 

As is well-known, A/"-fold SUSY is equivalent to weak quasi-solvability for one-dimensional 
Schrodinger operators and hence any of them admitting local solutions in closed form pos- 
sesses (at least) one type of A/"-fold SUSY. There have been four types so far found, and all 
the A/"-fold SUSY systems so far constructed are those having polynomial-type local solu- 
tions, or equivalently, those preserving a polynomial subspace. It should be noted, however, 
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that except for type A they constitute only subclass of the most general A/"-fold SUSY in 
each type. In the case of type A, it was already proved [2, 3] that it is equivalent to the 
sl(2) Lie-algebraic quasi-solvable models preserving a type A monomial space. On the other 
hand, all the other types of jV-fold SUSY systems were constructed under a particular re- 
striction F(q) = z'{q)/z{q) which guarantees that a linear function space preserved by the 
operator under consideration is of monomial or polynomial type. In other words, linear 
differential operators admitting an invariant subspace of non-polynomial type have not been 
systematically considered yet. 

To explore the area of operators with a non-polynomial invariant subspace, it is now 
evident that we must tackle a necessary and sufficient condition for other types of A/"-fold 
SUSY, or at least examine them under a less restrictive condition than F(q) = z'(q)/z(q). 
The present investigation on the necessary and sufficient condition for type B 3-fold SUSY 
would definitely provide us the first clue and foothold toward that research direction. 

The notion of shape invariance introduced in Ref. [18] is a well-known tool for investigat- 
ing and constructing an exactly solvable Schrodinger operator. Precisely speaking, it is a 
sufficient condition for solvability but is neither necessary nor sufficient for exact solvability. 
Later in Ref. [19], it was generalized to the notion of two- and multi-step shape invariance, 
keeping the latter aspect, namely, sufficiency for solvability, intact. Recently, it was shown 
[14] that any two-step shape-invariant model can be systematically investigated and con- 
structed in the framework of A/"-fold SUSY as a particular case of 2-fold SUSY with an 
intermediate Hamiltonian developed in Ref. [13]. Now that we appreciate the most general 
3-fold SUSY which is quasi-solvable in the strong sense, we are in a position to follow the 
latter program to uncover the area of three-step shape invariance. In the case of type A 3- 
fold SUSY, it was already shown [20] that the possible three-step shape-invariant potentials 
are all reducible, that is, they have ordinal shape invariance as well. Hence, we would be 
able to see whether there exists an irreducible three-step shape-invariant potential by using 
the results here. 

Appendix A: Quasi-solvable operators in type C 3-fold SUSY 

In Ref. [15], the set of quasi-solvable operators up to second order preserving the type 
C monomial subspace (4.1) was first investigated and four were found for an arbitrary 
Afi,N~2 £ N, cf., the general case (4-a) there. In our present case of M\ = 2 and A2 = 1, 
they reduce to (cf., (3.7a)-(3.7d) in Ref. [4]) 



In addition to the above, there exist two extra second-order operators preserving the three- 
dimensional type C monomial space (4.2). This is because three-dimensional monomial 
subspaces are special and it can also be regarded as the special case A (a) in Ref. [15]. The 
two additional operators are given by 




(Al) 




(A2) 
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and were not considered in Ref. [4]. 

The other set of quasi-solvable operators related with the above via 3-fold SUSY preserves 
the linear space V3" in (4.8) which is equivalent up to a similarity transformation to another 

type C monomial subspace (1, z, z 1_A ). Hence, it is evident that they are given by K^ c \\) = 
zjf\l - X)z- 1 and A^ C) (A) = z j\f' '(1 - A) z' 1 = +,0,-), cf., (3.29) in Ref. [4]. 
Explicitly, their forms read as 

*<<»(*) = 4-1, *S?(A) = *-'(4-l)(4-2 + A), 
-(c) /n _/ d n\f„ d ^(C)/- ( d A / d 



A-(A) =(,211,11, A-(A) = * ^ - 2 - 2 + A 



(A3) 



'(C)/ _ / d -\ / d 



K${\) = z x ~ x (z^--2)(z^--2 + X 



The last two operators were not considered in Ref. [4] . 



Appendix B: Quasi-solvable operators in type B 3- fold SUSY 

The set of quasi-solvable operators up to second order preserving the type B monomial 
subspace (4.12) was also investigated in Ref. [15] and six were found for an arbitrary Af EN. 
In the case of Af = 3, they reduce to (cf., (3.11a)-(3.11f) in Ref. [5]) 

t(b) _ d t(b) _ d 2 ( B ) _ d / d \ ( B ) _ 2 d 2 

<Jc\ — Z , , J — , _ , Jr,_ — , I Z , O I , J nr. — Z 



" 'dz' (I-'' d: \* dz J' "" " d: 2 ' . pn 



f " 1 d^ ; \ dz v ++ " d: V"<1- 

In addition to the above, there exists an extra second-order operator preserving the three- 
dimensional type B monomial space (cf., the special case A (b) in Ref. [15]): 

^W(4-3)(4-!). (B2) 

The latter operator was not considered in Ref. [5] . 

The other set of quasi-solvable operators related with the above via 3-fold SUSY, which 
preserves the linear space in (4.19), reads as (cf., (3.16a)-(3.16f) in Ref. [5]) 

(B3) 

K™ = z(z±-2\ (z± + l\, K?l = z*[z--2\[z- 



+0 1 <k \ dz r ++ ~ Vdz -jvdz 
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The space is equivalent to a three-dimensional monomial subspace (l,z 2 ,z 3 ) up to a 
similarity transformation. Hence, the latter monomial space is preserved by z K z and 
z z~ l = +, 0, — ). A set of quasi-solvable operators of up to second order preserving 
it was also examined in Ref. [15] as the special case A (c) with m = 2 and six were presented. 
Transformed back to the space V 3 + in (4.19), they correspond to K { *\ K^J, K { *\ K^, 
K+l, and another one given by 

^('aHW- 1 )- (B4) 

Hence, there have been in total seven linearly independent ones so far found in the litera- 
ture [5, 15]. 



Appendix C: Quasi-solvable operators in type A 3- fold SUSY 

The set of quasi-solvable operators up to second order preserving the type A monomial 
subspace (4.22) is exhausted by the first-order differential operator representation of sl(2) 
Lie algebra and their quadratic form, which was first reported in Ref. [17] (cf., (4.13) and 
(4.17)-(4.21) in Ref. [2]). 

J (A) - — J (A) - z— J (A) - z (z— - 
J ~ dz' J ° ~ dz' J+ - \ Z dz 

J (A) -— J^- Z f- j^- z ^ (CI) 
J ~ d^ 2 ' J °" d^ 2 ' J °° ~ Z dz 2 ' 

Jt$ = (z± - 2] , J&) = , 2 (z$- - 2] (z$- - 1 



+0 ~ dz V dz J ' ++ " \~dz 'J \'dz 

The other set of quasi-solvable operators related with the above via 3-fold SUSY preserves 
the same type A monomial subspace (4.22) and thus is the same, as was shown in Ref. [3] 
for arbitrary JVgN. Hence, we have = jf^ and = (i = +, 0, — ). 



Appendix D: Quasi-solvable operators in type X2 3- fold SUSY 



In Ref. [8], second-order quasi-solvable operators preserving a type X 2 exceptional poly- 
nomial subspace were first studied and four linearly independent ones were found for each 
subspace. In the case of three-dimensional subspace V3" [u; a] = Vp°V;a], namely, (5.1) 
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with (5.15), they reduce to (cf., (2.12), (2.15), (2.18), and (2.20) in Ref. [8]) 

ux->), \ d 2 , „. d 4(a — l)(u + a) ( d 

Ji ^^-(u-a + 3)- + f ^ a) 

4 X2 \a) = [u 2 + (a + 2)(a- 1)] -A - [u 2 + 4u + (a - l)(3a + 2)] A + 4« 

att + a 2 - 1 / d 

d 2 d 
A X2 \a) = (w + 2a + 2)u 2 — + [(a - 5)w 2 + (3a 2 + a - 8)u + 4(a + 4) (a - 1)1 — 

dir dw 

(a 2 + 3a-8)w + (a + 4)(a- l)a / d 

- 4 a - 2)u - 4 a - 1 y ' ^ '— — - 1 

/(u;a) \ dM 

d 2 

2(a + l)4 Xa) (a) = [2(a + l)u + 3a 2 + 7a + 6] w 3 — - [l2(a + l)u 3 

+ (3a 3 + 8a 2 + 15a + 22)w 2 + (a - l)(7a 3 + 27a 2 + 10a - lQ)u 

+ 2(a - l)(a 4 + 8a 3 + 29a 2 - 6a - 40)] A 



+ 24(a + l)w 2 + 4(a - l)(3a 2 + 2a - A)u 

0)w+ (a 3 

f(u;a) \du 



x9 (a 3 + 9a 2 - 22a - 40)w + (a 3 + 9a 2 - 6a - 20)a / d 
+ 4 (« - 1) 777^ XT " 1 



(Dl) 



The other set of quasi-solvable operators related with the above via 3-fold SUSY, which 
preserves the other X 2 exceptional polynomial subspace V^ X2b \u; a] defined in (5.20) with 
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(5.21), reads as (cf., (3.38), (3.41), (3.44), and (3.46) in Ref. [8]) 

d 2 d 4 d 
K[ X2 \a) = u-—r + (u - a - 3)- — h -77 r (a - l)(u + a)- — h a(u + a - 1) 

d-u 2 d-u /(-u; a) l du 



Au 



K { 2 X2 \a) = [u 2 + (a - l)(a - 4)] + [u 2 - 6u + (a - l)(3a - 8)] A - 

" 8 f*faa) { [( ° " 3)U + (« - - 2)]^ + (« - 2)w + « 2 - 3« + 4} , 

^( X2 )(a) = ( u + 2a - A)u 2 — - Ua + 3)u 2 + (3a 2 -5a- A)u - 4(a - l)(a - 2)1 — 

du 2 du 

+ Aau - ^ — ^ { [(a 2 - 3a + 4)it + a(ct - l)(a - 2)1 

/(■u; a) |^ L J d-u 

+ a[(a — 2)u + a(a — 3)]|, 

d 2 

2(a - 2)fsrf 2) (a) = [2(a - 2)u + 3a 2 - 11a + 12] m 3 — - [l2(a - 2)u 3 

- (3a 3 - 36a 2 + 97a - 84)m 2 - (a - l)(7a 3 - 49a 2 + 112a - 80)w 

d 

- 2(a - l)(a 4 - 11a 3 + 32a 2 - 36a + 16)1 — + 24(a - 2)m 2 

d-u 

- 4(3a 3 - 27a 2 + 64a - 48)w + - ^ { (a - 1) \{a 3 - 3a 2 + 8a - 16)u 

f{u;a) { 

d 

+ a(a 3 - 3a 2 + 8a - 16)1 — + a [(a 3 - 3a 2 + 8a - 16)u 

du 

+ a(a-l)(a 2 -3a + 4)]|, 



(D2) 
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